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Abstract 

i-^^l We prove a generalization of the quantum Markovian equation for observables. 

In this generalized equation, we use superoperators that are fractional powers of 
^ \ completely dissipative superoperators. We prove that the suggested superopera- 

^ \ tors are infinitesimal generators of completely positive semigroups and describe the 

properties of this semigroup. We solve the proposed fractional quantum Markovian 
equation for the harmonic oscillator with linear friction. A fractional power of the 
\ Markovian superoperator can be considered a parameter describing a measure of 

\^ , "screening" of the environment of the quantum system: the environmental influ- 

. ence on the system is absent for a = 0, the environment completely influences the 

Q>^ . system for a = 1, and we have a powerlike environmental influence for < a < 1. 

o 

O 



1 Introduction 

Fractional calculus appeared in 1695, when Leibniz described the derivative of order 
a = 1/2 [Il[2l[3]. Derivatives and integrals of noninteger order were studied by Leibniz, 
Liouville, Grunwald, Letnikov and Riemann. Many books have now been written about 
fractional calculus and fractional differential equations IH [21 HI [5l [6l [TJ [8] . Derivatives and 
integrals of noninteger order and fractional integro-differential equations have found many 
appUcations in recent studies in physics (see, e.g., [H [IDl [HI [12] and [T5 | [TH [TS| [TB]). 

In quantum mechanics, observables are given by self-adjoint operators. The dynamical 
description of a quantum system is given by superoperators. A superoperator is a map 
that assigns one operator some other operator. 
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The motion of a system is naturally described in terms of the infinitesimal change of 
the system. The equation for a quantum observable is called the Heisenberg equation. For 
Hamiltonian quantum systems, the infinitesimal superoperator is defined by some form 
of derivation. A derivation is a linear map L that satisfies the Leibnitz rule Ci^AS) = 
{CA)B + A{CB) for any operators A and B. A fractional derivative can be defined as 
the fractional power of the derivative (see, e.g., |l7j). It is known that the infinitesimal 
generator C = {l/ih)[H, . ], which is used for Hamiltonian systems, is a derivation of 
quantum observables. In [18], we regarded a fractional power of the derivative operator 
C = {1/ih) [H, . ] as a fractional derivative on a set of observables. As a result, we obtain a 
fractional generalization of the Heisenberg equation, which allows generalizing the notion 
of Hamiltonian quantum systems. We note that a fractional generalization of classical 
Hamiltonian systems was suggested in [19] (also see pO]). In the general case, quantum 
systems are non-Hamiltonian and C is not a derivation operator. For a wide class of 
quantum systems, the infinitesimal generator C is completely dissipative [2ll [221 ISSl l2l| • 
Therefore, it is interesting to consider a fractional generalization of the equation of motion 
for non-Hamiltonian quantum systems using a fractional power of a completely dissipative 
superoperator. 

The most general change of state of a non-Hamiltonian quantum system is a quan- 
tum operation [25l [26l EZl EH ESI EO]. A quantum operation for a quantum system can 
be described starting from a unitary evolution of some closed Hamiltonian system if the 
quantum system is a part of the closed system [321 133]. But situations can arise where it 
is difficult or impossible to find a Hamiltonian system that includes the given quantum 
system. As a result, the theory of non-Hamiltonian quantum systems can be consid- 
ered a fundamental generalization of the quantum mechanics of Hamiltonian systems 
[2T1 [221 [23l [23] . The quantum operations that describe the dynamics of non-Hamiltonian 
systems can be regarded as real completely positive trace-preserving superoperators on 
some operator space. These superoperators form a completely positive semigroup. The 
infinitesimal generator of this semigroup is completely dissipative. The problem of non- 
Hamiltonian dynamics is to obtain an explicit form for the infinitesimal generator, which 
is in turn connected with the problem of determining the most general explicit form of this 
superoperator. This problem was investigated in [331 ISSl EH] . Here, we consider super- 
operators that are fractional powers of completely dissipative superoperators. We prove 
that the suggested superoperators are infinitesimal generators of completely positive semi- 
groups. The quantum Markovian equations with a completely dissipative superoperator 
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are the most general form of the Markovian master equation describing the nonunitary 
evolution of a density operator that is trace preserving and completely positive. We con- 
sider a fractional generalization of the quantum Markovian equation, which is solved for 
the harmonic oscillator with friction. We can assume that other solutions and properties 
described in [371 EHl EHl llQl llll US SHI iU Il5| can also be considered for fractional gen- 
eralizations of the quantum Markovian equation and the Gorini-Kossakowski-Sudarshan 
equation |35] . 

A fractional power of infinitesimal generator can be considered a parameter describing 
a measure of "screening" of the environment. Using the interaction representation of the 
quantum Markovian equation, we consider a fractional power a of the non-Hamiltonian 
part of the infinitesimal generator. We obtain the Heisenberg equation for Hamiltonian 
systems in the limit as a — 0. In the case a = 1, we have the usual quantum Markovian 
equation. As a result, we can distinguish the following cases: (1) absence of the environ- 
mental influence (a = 0), (2) complete environmental influence (a = 1), and (3) powerlike 
screening of the environmental influence (0 < a < 1). The physical interpretation of the 
fractional quantum Markovian equation can be connected with an existence of a powerlike 
"screening" of the environmental influence. Quantum computations by quantum opera- 
tions with mixed states (see, e.g., [30]) can be controlled by this parameter. We assume 
that there exist stationary states of open quantum systems [46l [37l |47l HH |49l HQI |45] 
that depend on the fractional parameter. We note that it is possible to consider quantum 
dynamics with a low fractionality by a generalization of the method proposed in [50] (also 
see [5I1E2]). 

In Section 2, we briefly review of superoperators on an operator Hilbert space and 
quantum operations and introduce the notation. In Section 3, we consider the fractional 
power of a superoperator. In Section 4, we suggest a fractional generalization of the 
quantum Markovian equation. In Section 5, we describe the properties of the fractional 
semigroup. In Sections 6 and 7, we solve the fractional equations for the quantum har- 
monic oscillator with and without friction. 

2 Superoperator and quantum operations 

Quantum theories essentially consist of two structures: a kinematic structure describ- 
ing the initial states and observables of the system, and a dynamical structure describing 
the change of these states and observables with time. In quantum mechanics, the states 
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and observables can be given by operators. The dynamical description of the quantum 
system is given by a superoperator, which is a map from a set of operators into itself. 

Let M. be an operator space. We let Ai* denote the space dual to TW. Hence, A^* is 
the set of all linear functionals on Ai. The classic denotations for an element of Ai are 
\B) and B. The symbols {A\ and u denote the elements of Ai*. By the Riesz-Frechet 
theorem, any linear continuous functional uj on an operator Hilbert space At has the form 
u!{B) = {A\B) for all B G Ai, where 1^4) is an element in Ai. Therefore, the element 
A can be considered not only an element \A) of Ai, but also an element {A\ of the dual 
space At*. The symbol {A\B) for a value of the functional {A\ on the operator \B) is the 
graphic combination of the symbols {A\ and \B). 

Definition 1. A linear superoperator is a map C from an operator space Ai into itself 
such that the relation 

C{aA + bB) = aC{A) + bC{B) 
is satisfied for all A,Be D{C) C Ai, where D{C) is the domain of C and a, 6 G C. 

A superoperator £ assigns each operator A G D{C) the operator C{A). 

Definition 2. Let C be a superoperator on Ai. An adjoint superoperator of C is a 
superoperator A = C on Ai* such that 

{A{A)\B) = {A\C{B)) (1) 

for all B G D{C) C Ai and A E D{A) C M*. 

Let A4 be an operator Hilbert space and £ be a superoperator on Ai. Then {A\B) = 
Tr[A^B], and equation ([T]) becomes 

Tr[{A{A)yB] = Tr[A^C{B)]. 

If Ai is an operator Hilbert space, then by the Riesz-Frechet theorem. At and Ai* are 
isomorphic, and we can define the self-adjoint superoperators. 

Definition 3. A self-adjoint superoperator is a superoperator C on a Hilbert operator 
space M such that iC{A)\B) = iA\C{B)) for all A, B e D{C) C Ai and D{C) = D{C). 
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Let be a normed operator space. The superoperator L is said to be called bounded 
if ||£(y4)||_A^ < c||y4||;vj for some constant c and all A G M.. The value 

\\L{A)U 
||£|| = sup — —- 

is called the norm of the superoperator L. If is a normed space and £ is a bounded 
superoperator, then ||£|| = ||£||. 

In quantum theory, the class of real superoperators is the most important. 

Definition 4. Let Ai be an operator space and be an adjoint operator of A E ^A. A 
real superoperator is a superoperator C on M. such that 

[C{A)]^ = C{A^) 

for all A e D{C) C M and A^ G D{C). 

If £ is a real superoperator, then A = £ is real. If £ is a real superoperator and A 
is a self-adjoint operator = A E D{C), then the operator B = C{A) is self-adjoint. 
Then superoperators from a set of quantum observables j\4 into itself should be real. All 
possible dynamics of quantum systems must be described by a set of real superoperators. 

Definition 5. A nonnegative superoperator is a map C from M. into Ai, such that 
> for all A? = A'' A G D{C) C A4. A positive superoperator is a map C from M. 
into itself such that C is nonnegative and C{A) = if and only if A = 0. 

Let denote an operator algebra. A left superoperator corresponding to A G Al is 
a superoperator La on Ai such that LaC = AC for all C G A^. We can think of La as 
meaning left multiplication by A. A right superoperator corresponding to A G Al is a 
superoperator Ra on Ai such that RaC = CA for all C G Ai. 

The most general state change of a quantum system is a called a quantum operation 
[251 f26\ fI7\ UEl EO] • A quantum operation is described by a superoperator S that is 
a map on a set of density operators. If p is a density operator, then S{p) should also be 
a density operator. Any density operator pt = p(t) is a self-adjoint (pj = pt), positive 
[pt > 0) operator with unit trace {Trpt = 1). Therefore, for a superoperator S to be the 
quantum operation, the following conditions must be satisfied: 
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1. The superoperator ^ is a real superoperator, i.e. (^S{A)^ = S{A^) for all A. The 

real superoperator S maps the self-adjoint operator p to the self-adjoint operator 
S{p): {S{p)y=S{p). 

2. The superoperator ^ is a positive superoperator, i.e. S maps positive operators to 
positive operators: S{A'^) > for all A ^ or £{p) > 0. 

3. The superoperator ^ is a trace-preserving map, i.e. {I\£\p) = = 1 or 
8^1) = I. 

Moreover, we assume that the superoperator S is not only positive but also completely 
positive [53]. The superoperator £^ is a completely positive map from an operator space 
M into itself if 

n n 
k=l 1=1 

for all operators Ak, Bk E Ai and any integer n. 

Let the superoperator £^ be a convex linear map on the set of density operators, i.e. 

s s 

where < < 1 for all s, and Y2s = 1- Any convex linear map of density operators 
can be uniquely extended to a linear map on self-adjoint operators. We note that any 
linear completely positive superoperator can be represented by 

m m 

£ = Y,La,RaI-- S{p) = Y.A,pAl 

k=l k=l 

If this superoperator is trace-preserving, then 



Y,AlA, = I. 



k=l 

Because all processes occur in time, it is natural to consider quantum operations 
Sit^to) that depend on time. Let the linear superoperators Sit^to) form a completely 
positive quantum semigroup [M] such that 

j^8{t,to)=kt8{tM). (2) 



where Aj is an infinitesimal generator of tfie semigroup [Ml EH El] • Tfie evolution of a 
density operator p is described by 

£{t,t^)p{t^)= pit). 

We consider quantum operations S{t,tQ) witli an infinitesimal generator A such that the 
adjoint superoperator C is completely dissipative, i.e. 

C{AkAi) - C{Ak)Ai - AkC{Ai) > 0, 

for all Ai,...,An G D{C) such that AkAi G D{C). The superoperator C describes the dy- 
namics of observables of a non-Hamiltonian quantum system. The completely dissipative 
superoperators are infinitesimal generators of completely positive semigroups t > 0} 
that are adjoint to {^^1 t > 0}, where St = S(t, 0). 



3 Fractional power of a superoperator 

Let £ be a closed linear superoperator with an everywhere dense domain D{C) and a 
resolvent R{z, C) on the negative semiaxis, and satisfy the condition 



We note that 



The superoperator 



\\R{-z,C)\\<M/z, (z>0,M>0). (3) 
R{-z,C) = {zLi + C)-\ 



sinvra r dzz''-'R{-z,C)C 

71" Jo 



(4) 



is defined on D{C) for < a < 1 and is called a fractional power of the superoperator C 
[55l [56] . We note that the superoperator allows a closure. If a closed superoperator C 
satisfies condition then C'C^ = /:"+^ for a, /3 > 0, and a + (3 < 1. 

Let £ be a closed generating superoperator of the semigroup {$^1 t > 0}. Then the 
fractional power of C is given by 

r(-a) Jo 

which is called the Balakrishnan formula. 
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The resolvent for the superoperator can be found by the equation 

R{-z,C") = {zLi + C")-^ = 
^snvna r ^ — 

TT Jo -2 + COS ttq; + 

called Kato 's formula. It follows from this formula that the inequality 

\\R{-z,C")\\<M/z, {z>0), 

is satisfied with the constant M in inequality ([3]) for the superoperator C. It follows from 
the inequality 

\\zR{-z,C)\\ = \\z{zLj + Cy'^W < M 

for all ^ > that the superoperator z{zLj + C)^^ is uniformly bounded in every sector of 
the complex plane given by the relation | argz| < for not greater than some number 
vr — -0, {0 < ip < Tc). Then the superoperator zR{—z,C°') is uniformly bounded in every 
sector of the complex plane such that | arg2;| < (p for < vr — aip. 

Let £ be a closed generating superoperator of the semigroup {^t\ t > 0}. Then the 
superoperators 

POO 

$1"^ = / dsUt,s)^s, (t >0), (5) 
Jo 

form a semigroup such that is an infinitesimal generator of . Equation ([5]) is called 
the Bochner-Phillips formula. 

In equation ([5]), we use the function 

fait, = ^ / exp(sz -tz'^), (6) 

■^'^'^ J a— too 

where a, t > 0, s > 0, and < a < 1. The branch of is chosen such that Re{z°') > for 
Re{z) > 0. This branch is a one- valued function in the z plane cut along the negative real 
axis. This integral obviously converges by virtue of the factor exp(— t^;"). The function 
fait, s) has the following properties: 

1. For all s > 0, the function fait, s) is nonnegative: fait, s) > 0. 

2. We have the identity 



noo 

/ dsfait,s) = 1. 
^0 
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3. For t>0 and a; > 0, 

poo 

-tx° 



/ cise-'^7„(t,s) = e" 
Jo 



4. Passing from the integration contour in ([6]) to the contour crossing of the two rays 
r exp(— i6') and r exp(+z6'), where r G (0, cxd), and 7r/2 < 6* < tt, we obtain 



faiths) 

vr 



1 /""^ 

— / cir exp(sr cos 6' — tr" cos(q;6'))- 
Jo 

■ sin(sr sin 6^ — sin(a0) + 6'). (7) 
5. If a = 1/2, then 6 = it, and 

Mt,s) = -| dre-sin(tv^) = ^^3/^, 

which is a corollary of equation ([7]). 



4 Fractional quantum Markovian equation 

The motion of a systems is naturally described in terms of the infinitesimal change. 
This change can be described by an infinitesimal generator. One problem of the non- 
Hamiltonian dynamics is to obtain an explicit form of the infinitesimal generator. For 
this, it is necessary to find the most general explicit form of this superoperator. The prob- 
lem was investigated in [3ll [351 ESj for completely dissipative superoperators. Lindblad 
showed that there exists a one-to-one correspondence between the completely positive 
norm-continuous semigroups and completely dissipative generating superoperators [36] . 
Lindblad's structural theorem gives the most general form of a completely dissipative 
superoperator. 

Theorem 1. A generating superoperator Cy of a completely positive unity-preserving 
semigroup {$i = exp(— t£y)| t > 0} on an operator space M. can he represented in the 
form 

1 ^00 

-Cy{A) = --^[H,A] + + iVk,^]Vk), (8) 

k=l 

where H, Vk, Efc V^, V^Vk E M. 
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We note that the form of Ly is not uniquely fixed by ([H]). Indeed, formula (IHl) preserves 
its form under the changes 

^ oo 

Vk ^ H + a,/, E ^ ^ + 

fc=i 

where are arbitrary complex numbers. 

Using At = where = exp(— t£y), we obtain the equation 

7 1 -.00 

fc=l 

where £y is defined by ([8]). This is called the quantum Markovian equation for the 
observable A. 

The Lindblad theorem gives an explicit form of the equations of motion if the following 
restrictions are satisfied (here Ay is adjoint to £y): 

1. £y and Ay are bounded superoperators and 

2. £y and Ay are completely dissipative superoperators. 

Davies extended the Lindblad result to a class of quantum dynamical semigroups with 
unbounded generating superoperators [57] . 

We consider quantum Markovian equation (Q for an observable At. We rewrite this 
equation in the form 

^At = -CvAt, (10) 
at 

where £y denotes the Markovian superoperator 

00 

fc=i 

Here, we use the superoperators of left multiplication Ly and right multiplication Ry 
determined by the relations LyA = VA and RyA = AV. The superoperator Ljj is a left 
Lie multiplication by A such that 

LJ,A=^[H,A]. (12) 

If all operators Vk are equal to zero, then Cy = Ljj, and equations (fTOl) and ffTTl) give the 
Heisenberg equations for a Hamiltonian system. In the general case, the quantum system 
is non-Hamiltonian [2l] . 
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We obtain a fractional generalization of the quantum Markovian equation. For this, 
we define a fractional power for the Markovian superoperator Ly in the form 

Sin Trrv z**^ 

-[^CyY = - -\ dzz''-^R{~z,Ly)Ly (0<a<l). (13) 

The superoperator i^Cy)'^ is called a. fractional power of the Markovian superoperator. We 
note that {Cy)"'{Cv)^ = (Cy)"''^^ for a, P > 0, and a + (3 < 1. As a result, we obtain the 
equation 

j^A, = -{CyrA, (14) 

where t, H/h and Vk/\/h are dimensionless variables. We call this is the fractional quan- 
tum Markovian equation. 

If Vk = 0, then equation f|T4l) gives the fractional Heisenberg equation [18] of the form 

j^A, = -{L-^TA. (15) 

The superoperator (Ljj)"' is a fractional power of the left Lie superoperator fll2p . We note 
that this equation cannot be represented in the form 

d _ i 

— At = -L^^^^At = -[Hnew,At\ 

with some operator Hmw Therefore, quantum systems described by (ITSll are not Hamil- 
tonian systems. These systems are called the fractional Hamiltonian quantum systems 
(FHQS). Usual Hamiltonian quantum systems can be considered a special case of FHQS. 
We note that a fractional generalization of classical Hamiltonian systems was suggested 

in ^m\. 

Using the operators 

Au{t) = U{t)AtU\t), Wk{t) = U{t)VkU\t), 
where f/(t) = exp{(l/-i/i)/7}, we can write the quantum Markovian equation in the form 

^Auit) = -CwAu{t). (16) 
at 

The superoperator 

. oo 

k=l 
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describes the non-Hamiltonian part of the evolution. Equation (JT6l) is the quantum 
Markovian equation in the interaction representation. The fractional generalization of 
this equation is 

j^Auit) = -{CwTAuit). (18) 

Equation (fTS!) is the fractional quantum Markovian equation in the interaction represen- 
tation. The parameter a can be considered a measure of the influence of the environment. 
For a = 1, we have quantum Markovian equation ( |T6l) . In the limit as a — 0, we obtain 
the Heisenberg equation for the quantum observable At of a Hamiltonian system. As a 
result, we can consider the physical interpretation of equations with a fractional power of 
the Markovian superoperator an influence of the environment. The following cases can 
be considered in quantum theory: (1) absence of the environmental influence {a = 0), 
(2) complete environmental influence (a = 1), and (3) powerlike screening of the environ- 
mental influence (0 < a < 1). The physical interpretation of fractional equation (ITS]) can 
be connected with an existence of a powerlike screening of the environmental influence on 
the system. 



5 Fractional semigroup 

If we consider the Cauchy problem for equation (ITOl) with the initial condition given 
at the time t = by Aq, then its solution can be written in the form At = ^tAo. The 
one-parameter superoperators t > have the properties 

= $t+„ {t,s>0), ^Q = Li. 

As a result, the superoperators form a semigroup, and the superoperator Cy is a 
generating superoperator of the semigroup t > 0}. 

We consider the Cauchy problem for fractional quantum Markovian equation f[T^ with 
the initial condition given by Aq. Then its solution can be presented in the form 

At{a) = <I>1")Ao, 

where the superoperators t > 0, form a semigroup, which we call the fractional 

semigroup. The superoperator — (£y)" is a generating superoperator of the semigroup 
{$^"^1 t > 0}. We consider some properties of the fractional semigroups | t > 0}. 
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The superoperators $J can be constructed in terms of $f by Bochner-Phillips formula 
([5]), where fa{t,s) is defined in ([6]). If At is a solution of quantum Markovian equation 
(ITU]) , then formula (jS]) gives the solution 

POO 

At{a) = / dsf^{t,s)As, {t > 0) 
Jo 

of fractional quantum Markovian equation ([H]). 

A linear superoperator ^j""* is completely positive if 

J2B^<^i"\AlA,)B,>0 

for any Ai,Bi G The following theorem states that the fractional semigroup {^["^ \ t > 
0} is completely positive. 

Theorem 2. t > 0} is a completely positive semigroup of superoperator $t on Ai, 

then the fractional superoperators form a completely positive semigroup | t > 0} . 

Proof. Bochner-Phillips formula gives 

dsf^it,s)Y,B^^siAlA,)B, 

for t > 0. Using 

J2Bi<^siAlA,)B,>0, fa{t,s)>0 (s>0), 

we obtain 

id 

Corollary. If ^t, t > 0, is a nonnegative one-parameter superoperator, i.e., > 

for A>Q, then the superoperator (^'f^ is nonn egative, i.e., $J"^(A) > /or A > 0. 

Using the Bochner-Phillips formula and the property fa(t,s) > 0, s > 0, we can 
easily prove that the superoperator is nonnegative, if t > is a nonnegative one- 
parameter superoperator. This corollary can also be proved by using Bi = I, Ai = A, 
and Ai = Bi = {i = 2, ...) in the proof of the theorem. 
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In quantum theory, the class of real superoperators is the most important. Let 
G Ai* be adjoint to A & Ai. A real superoperator is a superoperator $j on A^, 
such that {^tAy = ^t{A'^) for all A G -D($i) C M. A quantum observable is a self- 
adjoint operator. If $t is a real superoperator and A is a self-adjoint operator, A'^ = A, 
then the operator At = ^tA is self-adjoint, i.e., {^tAY = ^tA. Let A4 he a. set of quan- 
tum observables. Then superoperators on into must be real because quantum 
dynamics, i.e., temporal evolutions of quantum observables, must be described by real 
superoperators. 

Theorem 3. If^t is a real superoperator, then the superoperator ^["^ is also real. 
Proof. The Bochner-Phillips formula gives 

POO 

($^)A)t= / dsat,s){^sA)\ (t>0). 
Jo 

Using (jTj), we can easily see that fa(t,s) = fa(t,s) is a real- valued function. Then 
($,A)t = $iAt leads to (<I>^^A)t = <!>i"\A^) for all A e D{<^1°'^) C M. □ 

If $i is a superoperator on a Hilbert operator space Ai, then an adjoint superoperator 
of is a superoperator St on J^* such that 

{St{A)\B) = {A\MB)) (19) 

for all B G -D($t) C Al and some A G Al*. Using the Bochner-Phillips formula, we 
obtain the following theorem. 

Theorem 4. If St is an adjoint superoperator of^t, then the superoperator 

POO 

Si''^ = / dsf^it,s) Ss, (t>0), 
Jo 

is an adjoint superoperator of ^[""^ . 

Proof. Let St be adjoint to i.e. equation f|T9|) is satisfied. Then 

POO 

iSi''^A\B)= / dsf^{t,s){SsA\B) = 
Jo 
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/■oo 

/ dsUt,s)iA\^sB) = {A\^l''^B). □ 
Jo 



It is known that St is a real superoperator if $j is real. Analogously, if is a real 
superoperator, then Sj:°'^ is real. 

Let {St\t > 0} be a completely positive semigroup such that the density operator 
Pt = StPo is described by 

j^pt = -AvPt, (20) 

where Ay is adjoint to the Markovian superoperator Cy The superoperator Ay can be 
represented in the form 

1 oo 

AvPt = --^[H,Pt] + ;^ 5^(l4.p*^J - (ptV^Vk + V^VkPt) 

k=l 

We note that equation (!20|) with = gives the von Neumann equation 

^.Pt = ^[H,pt\- 
at in 

The semigroup {Sj: ''| t > 0} describes the evolution of the density operator ptiot) — £\: ^ po 
by the fractional equation 

^pt(a) = -{kvTpt{a). 

This is the fractional quantum Markovian equation for the density operator. For = 0, 
this equation gives 

which can be called the fractional von Neumann equation. 



6 Fractional equation for the harmonic oscillator 

We consider a quantum harmonic oscillator such that 

H = ^P'+'^Q\ V, = 0, (21) 
2m 2 

where t and P are dimensionless variables. Then equation (fUll (also see ( ITSl) ) describes 
a harmonic oscillator. For A = Q and A = P, equation (HM for a = 1 gives 

^Qt = -Pt, ^Pt = -muj^Qt. 
at m at 
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The well-known solutions of these equations are 

Qt = Qo cos(ut) + ^— Po sin(co't), 

muj 

Pt = Pq cos{ujt) — muQQ sin(to't). (22) 

Using these solutions and the Bochner-Phillips formula, we can obtain solutions of the 
fractional equations 

j^Qt = -{L7,rQu JtP^ = -^LHTPt. (23) 
where H is given by fl^Tl) . The solutions of fractional equations fl2^ have the forms 

Jo 

POO 

Pi(«) = <l>S")Po = / dsUt,s)Ps. (24) 
Jo 

Substituting (122!) in gives [18J the equations 

Qt = QoCait) + —PoS^it), Pt = PoCa(t) - mcuQoSait), (25) 

where 

/•OO /"OO 

Cait) = / ds fa{t, s) cos{ujs), Sa{t) = / rfs /q (t, s) sin(ci;s) . 



Equations (l25l) describe solutions of fractional equations fl23|) for the quantum harmonic 
oscillator. For a = 1/2, we have 



These functions can be represented in terms of the Macdonald function (see Sec. 2.5.37.1 
in [58]), which is also called the modified Bessel function of the third kind. 
It is easy to obtain the expectations 

<Qt>= XoCo,{t) + —poSa{t), 

muj 

< Pt >= PoCa(t) - mUJXoSait), 
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and the dispersions 

Here, we use the coordinate representation and the pure state 

^{x) =< >= —L= exp{-^^^^ + '-Pox}. (26) 
The expectation and dispersion are defined by usual 



7 Fractional quantum Markovian equation for the os- 
cillator with friction 

We consider the fractional quantum Markovian equation with 7^ 0. The basic 
assumption is that the general form of a bounded completely dissipative superoperator 
given by the quantum Markovian equation also holds for an unbounded completely dissi- 
pative superoperator Cy- Another condition imposed on the operators H and Vk is that 
they are functions of the operators Q and P such that the obtained model is exactly 
solvable ES] (also see SD])- We assume that Vk = Vk{Q,P) are the first-degree 
polynomials in Q and P, and that H = H{Q, P) is a second-degree polynomial in Q and 
P. These assumptions are analogous to those used in classical dynamics when friction 
forces proportional to the velocity are considered. Then Vk and H are given in the forms: 

H = —P^ + -—Q^ + ^{PQ + QP), Vk = akP + bkQ, (27) 
where and 6^, k = 1,2, are complex numbers. It is easy to obtain 

CvQ = -P + fiQ- \Q, 
m 

CyP = -mu^Q - fiP- \P, 

where 

n=2 n=2 

A = Im(^akb*k^ = -Im(^^a*kbky 

k=l k=l 
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Using the matrices 



m 



PI \ —muj^ — \ 

we write the quantum Markovian equation for At as 



at 



where CyAf = MAf. The solution of (EHD is 



OO CO ^ 

A, = $,Ao = V -^C^Ao = V ^MMo. 

n=0 n=0 

The matrix M can be represented in the form M = N~^FN, where F is a diagonal 
matrix. Let i/ be a complex parameter such that u"^ = fi"^ — u"^. Then we have 

-(A + z/) 
-{X-u) 



Taking 

n=0 ■ \n=0 

into account, we obtain the superoperator $f in the form 



cosh(z/t) + (/i/z^) sinh(z/t) (1/mz/) sinh(z/t) 

— (mc<j^/z/) sinh(z/i(;) cosh(z/t) — (/x/i/) sinh(i/t) 

As a result, we obtain 

Q^ = e"^*[cosh(z/t) + - sinh(z/)f:)]Qo + — e"^* sinh(z/t)Po, 

2 

Pt = —^'^^ sinh(z/t)Qo + e"^* [cosh(z/t) - ^ sinh(i/t)]Po- (29) 

The fractional quantum Markovian equations for Qt and Pt are 

= -(/:yrgt, AQ^ = -(/:v'rQ„ (so) 
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where t and Vk/Vh are dimensionless variables. The solutions of these fractional equations 
are given by the Bochner-Phillips formula, 



Qt{a) = <l>l"^Qo = / dsfc.it, s)Qs, (t > 0), 



POO 

Pt{a) = <l>f^Po = / dsf^it, s)Ps, (t > 0), 
Jo 



(31) 



where Qg and Pg are given by (l29l) and the function fa{t, s) is defined in IQ. Substituting 
(I29D in daH) gives 

1 



Q,{a) = [Ch^it) + ^SK{t)]Q, + — C/i„(t)Po, 



Ptia] 



-SK{t)Q, + [CK{t) - ^SK{t)]Po 



(32) 



where 



Cha{t) = / ds fa{t,s) e~^^ coshius), 
Jo 



Sha{t) = / dsfa{t,s)e '*sinh(z/s). 



For a = 1/2, we have 



t 



t 



COSh(t/g) ^2/4,.;,, 

"■^ „3/2 ^ ' 



sinh(;/g) t2/4,_;,. 



These functions can be represented in terms of the Macdonald function (see Sec. 2.4.17.2 
in [HH]) such that 

\v{t,\,-u) + V{t,X,u) 



where we use the notation 

Vit,X,u) = 



2v^ 

t 



V{t,X,-u)-V{t,X,u) 



2v^ 



4A 



^-1/2(2 



where Reif') > Re{u), Re{X) > 0, and Ka{z) is the Macdonald function P[2]. 

As a result, equations (132!) define a solution of the fractional quantum Markovian 
equation for the harmonic oscillator with friction. 
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8 Conclusion 



Quantum dynamics can be described by superoperators. A map assigning each operator 
exactly one operator is called a superoperator. It is natural to describe motion in terms 
of the infinitesimal change of a system. The equation of motion for a quantum observable 
is called the Heisenberg equation. For Hamiltonian quantum systems, the infinitesimal 
superoperator is some form of derivation. A linear map L satisfying the Leibnitz rule 
Li^AS) = {CA)B -\- A{CB) for all operators A and B is called derivation. It is known that 
the infinitesimal generator C = {l/ih)[H, . ], which is used for Hamiltonian systems, is a 
derivation of quantum observables. We can regard a fractional power of the derivative 
C = {l/ih)[H, . ] as a fractional derivative on a set of quantum observables [H]. As a 
result, we obtain a fractional generalization of the Heisenberg equation [18], which allows 
generalizing the notion of Hamiltonian quantum systems. In the general case, quantum 
systems are non-Hamiltonian and C is not a derivation. For a wide class of quantum 
systems, the infinitesimal generator C is completely dissipative [2T], [221 ESI El] ■ 

Here, we consider a fractional generalization of the equation of motion for non-Hamiltonian 
quantum systems using a fractional power of a completely dissipative superoperator. We 
suggested a generalization of the quantum Markovian equation for quantum observables. 
In this equation, we used a superoperator that is a fractional power of a completely dis- 
sipative superoperator. We proved that the suggested superoperator is an infinitesimal 
generator of a completely positive semigroup and described properties of this semigroup. 
We solved the proposed fractional quantum Markovian equation exactly for the harmonic 
oscillator with linear friction. A fractional power a of the quantum Markovian superopera- 
tor can be considered a parameter describing a measure of "screening" of the environment. 
We can separate the cases where a = 0, absence of the environmental influence; where 
a = 1, complete environmental influence; and where < a < 1, a powerlike environ- 
mental influence. A one-parameter description of a screening of the coupling between 
the quantum system and the environment is thus a physical interpretation of a fractional 
power of the quantum Markovian superoperator. 

We note that the quantum Markovian equation describes a coupling between a quan- 
tum system and an environment (see [32]). Another physical interpretation of a fractional 
power of the inflnitesimal generator is connected with Bochner-Phillips formula as fol- 
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lows. Using the properties 

POO 

/ Ut,s) = l, fait,s)>0 {for all s > 0), 
Jo 

we can assume that fa{t,s) is the density of a probabihty distribution. Then Bochner- 
Philhps formula can be considered a smoothing of the evolution $t with respect to 
the time s > 0. This smoothing can be considered a screening of the environment of the 
quantum system. 

The function fa{t, s) can be represented as the Levy distribution using a reparametriza- 
tion. We note that Levy distributions are solutions of fractional equations (see, e.g., 
[T3 l [59l [60l [6T] ) that describe anomalous diffusion. It is known that quantum Markovian 
equations are used to describe the Brownian motion of quantum systems [37j. Perhaps, 
the fractional generalization of quantum Markovian equations can be used to describe 
anomalous processes and random walks fl3[ [T5l [16] in quantum systems. 
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